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ABSTRACT 

We study the effects of gravitational lensing by galaxy clusters of the background of dusty star- 
forming galaxies (DSFGs) and the Cosmic Microwave Background (CMB), and examine the implica- 
tions for Sunyaev-Zcl'dovich-based (SZ) galaxy cluster surveys. At the locations of galaxy clusters, 
gravitational lensing modifies the probability distribution of the background flux of the DSFGs as 
well as the CMB. We find that, in the case of a single- frequency 150 GHz survey, lensing of DSFGs 
leads to both a slight increase 10%) in detected cluster number counts (due to a ~ 50% increase 
in the variance of the DSFG background, and hence an increased Eddington bias), as well as to a rare 
(occurring in ~ 2% of clusters) "fiUing-in" of SZ cluster signals by bright strongly lensed background 
sources. Lensing of the CMB leads to a 55% reduction in CMB power at the location of massive 
galaxy clusters in a spatially-matched single-frequency filter, leading to a net decrease in detected 
cluster number counts. We find that the increase in DSFG power and decrease in CMB power due 
to lensing at cluster locations largely cancel, such that the net effect on cluster number counts for 
current SZ surveys is sub-dominant to Poisson errors. 

Subject headings: gravitational lensing — galaxies: luminosity function, mass function — galaxie clus- 
ters: abundances — methods: numerical 



1. INTRODUCTION 

Galaxy clusters are the largest collapsed objects in the 
universe. In hierarchical scenarios of structure formation 
their abundance strongly depends on the growth rate of 
structure over cosmic time, and measurements of their 
abundance as a function of mass and redshift can i mpose 
stringent cons traint s on cosmological p arameters (Wang 
fc Steinhard"t||1998[ [Holder et al.||200T| . 

'I'he majority oi baryonic mass m clusters is in the 
hot intra-cluster medium. In addition to emitting X- 
rays, this plasma interacts with Cosmic Microwave Back- 
ground (CMB) photons via inverse Compton sc attering, 
causing a local sp ectral distortion in the CMB ( Sunyaev 
&: Zeldovich|1972 ), known as the S unyae v- Z eldovich ( S Z ) 
ettect. I'his effect can be used to detect clusters in large 
area CMB surveys. The selection of clusters based on this 
method is particularly valuable since the surface bright- 
ness of the SZ effect is redshift-independent, and can 
therefore be used to explore a large range of redshifts. In 
addition the SZ signature depends stron gly on the total 
mass of the clusters wit h small scatter (Barbosa et al. 
1996[ Shaw et al. 2008), allowing surveys to construct 
nearly mass-limited catalogs 



Using data froni the South Pole Telescope (SPT) ( |Carl- 
strom et al. 2011), the first discover y of galaxy clusters 
via th e ti'L effect was reported in Staniszewski et al. 
(2009), followed by the first cosmological constraints 



from a cata log of SZ-detect ed clusters in Vanderlinde 
( |2010[ ). More recently [Marriage et al.| (|201ip re- 



et al. 



ported clusters discovered by the Atacama Cosmology 
Telescope (AC T), with cosmologic al results based on 



that catalog in|Sehgal et al. (2011), the Planck coUab 



oration released a sample oTSZ-se lect e d galaxy clusters 



(Planck Collaboration et al. 20111 and Reichardt et al 
(2012bp published a larger catalog from multi- frequency 
SFT 3ata. 

In SZ surveys, dusty star- forming galaxies (DSFGs), 
seen as an unresolved background of point sources, are 



a source of astronomical confusion which acts as a noise 
term in th e SZ-detectiqn of galaxy clusters; the n um- 
ber counts ([Vieira et al. || 2010[ ) and power spectrum ( Re- 
ichardt et al.||2012a|) ot these sources are well-studied. 



I'he other main source of confusion noise is the CMB, a 
Gaussian field with a n extremely well-characterized spa- 
tial power spectrum ( Keisler et al.|20lT Das et al.|2011 1 
on the scales of interest. 

At the location of galaxy clusters, due to the efficiency 
of clusters as gravitational lenses, these ba ckgrou nds can 
be strongly distorted, as emphasized by Blain (|l998). 
These distortions can significantly affect the probabil- 
ity distributions of these sky noise backgrounds, which 
can in turn affect the SZ detection statistics of such 
surveys and the resul ting cosmolo gical parameter con- 
straints. For example, Blain (1998) found that the con- 



fusion noise from point sources can be amplified by a 
factor of three within an aperture comparable to the in- 
ner region of the cluster where strong lensing (multiple 
images) can occur. In this paper we simulate gravita- 
tional lensing of background DSFGs and the CMB by 
galaxy clusters, measure the probability distributions of 
the lensed fields, and estimate the effects of this lensing 
on SZ cluster surveys. 

In section fj2] we develop an analytical analysis of the 
non-Gaussianity of the background confused galaxies and 
the CMB due to cluster lensing, focusing on understand- 
ing the underlying processes that lead to this effect. In 
section |3] we describe the ray-tracing simulations that 
are performed to measure the shape of the background 
probability distributions for realistic cluster mass profiles 
and backgrounds. In section f|4] we study the impact of 
this background lensing on SZ cluster surveys, specifi- 
cally a single-band SPT-like survey. Discussion and con- 
clusion are presented in section f|5] In what follows we 
assume a spatially flat universe described by the best-fit 
WMAP7 ( )Komatsu et al.|2011 ) cosmological parameters 
[h = 0.71, CT8 = 0.801, SIm = 2669). 
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2. ANALYTICAL ANALYSIS OF BACKGROUND 
FLUCTUATIONS 

2.1. Uniform Magnification 

The simplest example of background lensing is the case 
of uniform isotropic magnification of the CMB by a small 
factor of fio (close to 1) in each d i rectio n of a region 
of sky. As shown in 



(2012 1, such a uni- 



Bucher et al. 

form magnification leads to a shift in the power spectrum 
Sci/ci = ^oidlnci/dln£ + 2) for small /Xo; for larger mag- 
nifications and a power-law power spectrum ce oc the 
variance will be scaled by /io"''"- For constant power per 
In interval (q oc there is therefore no shift, but one 
could expect a substantial difference for either the damp- 
ing tail of the CMB, where the spectrum is substantially 
steeper than or for a Poisson distribution of sources 
where ci is constant. 

For a Poisson distribution of sources, the power would 
scale as '-^ /i^, while in the CMB damping tail, where 
the power- law can be steeper than the power could 
scale more steeply than ~ fi~^. These effects therefore 
go in opposite directions and are likely to be comparable 
in amplitude. 

2.2. Discrete sources 

An analytical model for the flux distribution inside an 
aperture due to point sources can be developed for sim- 
ple cases of lensing, as well as in the absence of lensing. 
This model is helpful in that it provides a test of the 
numerical implementation as well as providing some in- 
tuition for the processes that influence the shape of the 
final distributions. 

In this section we assume the simplest conditions: a 
top-hat filter, and only consider sources at a single red- 
shift of Zg — 2.0. The differential number counts dN/dS 
are assumed to follow the dist ribution predicted b y the 
Durham semi-analytic model fBaugh et al. 2005), but 
with all sources placed at a single z. 

2.2.1. Unlensed DSFG Background 

In the absence of lensing the number of point sources 
with flux between Si and Si + AS falling inside an aper- 
ture of angular radius 0ap is simply proportional to the 
number counts dN/dS{Si) and the aperture area ttB^ 



ap ■ 
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(1) 



To get the flux distribution wc flrst consider only 
sources with flux Si. The distribution of the number 
of these sources is given by a Poisson distribution with a 
mean oi < Ni >. The probability of observing k sources 
each of flux Si in a given area of sky is then 



P^ik) 



kl 



-<Ni> 



(2) 



Since each source has a flux of Si, in terms of the total 
aperture flux S = kSi due to sources with flux Si, we can 
write the probability distribution (for discrete values of 

s/s,y. 

< AJ. >(S/S.) 



The probability distribution of the addition of two such 
source populations (with different fluxes) is given by the 
convolution of the probabilities Pij{S) = {Pi * Pj){S) 
and if all i — 0,...N flux bins are included the total 
probability is 

PTotal{S) = {Po*Pl*.-*PN){S) (4) 

2.2.2. Lensed DSFG Background 

In the presence of lensing, different regions of sky un- 
dergo different magniflcation factors, slightly complicat- 
ing the analysis. For s implicity, here we consider an 



SIS halo (Kormann et al. 1994) lensing unresolved point 
sources before moving onto detailed numerical simula- 
tions using a more complex lens and source model. 

If the aperture is centered on the halo, we can break 
the aperture area into M circular rings with equal widths 
A6 such that 



M 



M 



(lap = XI =51 ^'^^J 



(5) 



where ^ap is the total area of the top-hat aperture and 
9j is the angular radius of each bin. 

For SIS halos with Einstein radius 6e the magniflcation 
at position 9j is given analytically by 



\ej\-eE 



(6) 



We can calculate the flux distribution inside such a cir- 
cular ring j. The lensed flux of each population i is now 
Si J — fijSi (where the hat refers to lensed quantities). 
Lensing also changes the observed solid angles, diluting 
the source number counts; objects inside this ring truly 
reside in a region of the sky with area flj = ^j/fJ-j in the 
absence of lensing. Therefore the number of sources in 
the j'th ring is 



< N,,, >-- 



dN 

Is 



{S.OAS, 



21x9, A9j 



(7) 



The distribution of the total lensed flux (S — kfijSi) 
coming from each source population is again controlled 
by the Poisson distribution with a mean of < Ni^j > (for 
discrete values of S/fijSi) 



(8) 



{S/ifi,S,))\ 

The distribution of total flux inside the j'th ring is 

Pj,TotiS) = (P,=lj * P^^2J * ... * P^=N.J){S) (9) 

However since all the light inside the aperture is com- 
bined, the distribution of the total flux inside the aper- 
ture is again a convolution of the probabilities inside each 
ring 

PTotal{S) = [fi=lj = l * ... * Pi=N,o = l] * 

... * [P^=l,J = M * ... * ^.=Wj = m] {S) (10) 

The combined flux is shown in Figure [T] compared with 
a numerical simulation of this same configuration, show- 
ing excellent agreement. The main effects of lensing are 
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Fig. 1. — Flux distribution inside a top-liat filter with ra- 
dius Rfiiter = 0.5' centered on a circular SIS lens halo with 
M=10"'^^ Mq from point sources at redshift Zs ~ 2.0. The Gaus- 
sian and the skewed curves show the unlensed and the lensed flux 
distributions respectively. The solid (red) curves show the analyt- 
ical predictions and the dashed (black) curves show the results of 
simulations. 

to shift the peak of the PDF to lower fluxes, to broaden 
the distribution and to populate a high flux tail. The 
peak of the distribution has moved to lower fluxes as 
there is a high probability of magnifying a spot on the 
sky that does not host a source; the distribution has 
broadened because there is now more shot noise because 
the effective area being sampled is smaller with lower 
source number density; and the high flux tail has been 
populated by strong lensing events where a source hap- 
pens to fall in a region of high magnification. Lensing 
leaves the mean of the distribution unaltered. 

This calculation can be generalized to include other 
filters by weighting each region's distribution appropri- 
ately. However since fi{9) is generally not analytic for 
more complicated lens profiles, a numerical ray-tracing 
approach is preferred. In addition the above treatment 
assumes point sources, while in reality the non-zero size 
of sources has a sub stantial effect on lensing statistics 
(Hezaveh et al. 2012). For these reasons, we now turn to 
numerical simulations based on ray-tracing. 

3. SIMULATIONS 



We use the ray-tracing code of Hezaveh & Holder 
(2011) to simulate realistic maps of the lensed back- 
ground sources. As a lens model, we assume that the 
halo density pro file in the outer regio ns follows a spheri- 
cal NFW profile ( [Navarro et al.|1997l ) and include an SIS 
profile at the halo center to mimic the enhanced baryon 
density at the centers of galaxy cluste rs. We assume 
a mas s-concentration relation from Ma ndelbaum et al.l 
((2008) and a central velocity dispersion of 250 km/ s for 
tne SIS component. 

We generate lensing defiection maps for multiple source 
redshifts (ranging from 0.7 to 4.5 for DSFGs, and one at 
z = 1100 for the CMB). For DSFGs we populate the 
source plane with circular galaxies of radius 1 kpc, with 



numb er counts following the Durham model ( jBaugh et al. 
20051). For each source fiux bin we include i\ sources of 



the given flux in the source plane map, where N is picked 
randomly from a Poisson distribution with a mean of 
(dN/dS/dA) 6S Arnap where SS is the width of the fiux 
bins and Amap is the total area of the source plane grid. 



For lensing of the CMB we make realizations of random 
CMB fields with p ower spectra calculated using CAMB 
( [Lewis et all2000[ ). 

We thus generate and lens separate maps of the back- 
ground DSFGs and the CMB. An example of a lensed 
DSFG map is shown in Figure [2| It is clear from this 
figure that lensing can amplify the flux of a few highly 
magnified DSFGs while diluting their source density. 



t V •■ • ' : -Uf- ' ■ . 
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Fig. 2. — The inner 5x5 arcmin of a simulated field of DSFGs 
for a lens of 1 X IO-'^^Mq at = 0.5. As it can be seen from 
the figure lensing of DSFGs can both magnify the sources and the 
empty regions of the sky. 



4. EFFECTS ON SZ CLUSTER SURVEYS 

The impact of lensing on SZ cluster surveys will depend 
on the details of the experiment. Multifrequency surveys 
can isolate the SZ signal spectrally, while an instrument 
with a large beam will be mainly sensitive to regions 
of the cluster with relatively low magnification. As a 
concrete example, we will consider a si ngle-band 150 GHz 
cluster survey modelled after that of jVanderlinde et al.[ 
(20101 to investigate the impact of lensing. 



4.1. Matched-filter 

SZ surveys have generally employed a matched-filter 
approach to identifying galaxy clusters. The raw mi- 
crowave sky maps are filtered to optimize detection of 
objects with morphologies similar to the SZ signatures 
expected from galaxy cluster s, through the applicati on of 
spatial mat ched-filte rs (Haehnelt fc Tegmark ^^ l996 | Her- 
ranz et al.[[2002b|a[ [Melin et al.||2006[ ). I'br the single^ 
frequency case, in the spatial l^burier domain the map is 
multiplied by 



V'(fca;, ky) 



B{k,,ky)S{\k\) 



B{kx, ky)'^Nastro{\k\) + N„oise{kxj ky) 



where ip is the matched-filter, B is the response of the 
instrument after timestream processing to signals on the 
sky, S is the assumed source template, and the noise 
power has been broken into astrophysical (Nastro) and 
noise (Nnoise) components. For the source template, we 
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assume AT cx (1 + 9"^ /Ol) ^, where the core radius 6c is 
as sumed here to be 0.5', co nsistent with the resuhs found 



Reichardt et al. 



( |2012bD . 

i'his filter application is a linear operation, so the effect 
on the DSFG background and the CMB can be consid- 
ered independently. We simulate many realizations of 
the DSFG field and the CMB (pairs of lensed and un- 
le nsed maps), apply a matc hed-filter similar to that used 



Vanderlinde et al. (2010) to them, and record the val- 



ues of the filtered maps at the center of the cluster 
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Fig. 3. — The ratios of the standard deviations of the lensed 
and unlensed fields at the center of the cluster after convolution 
with the matched-filter. For each mass and redshift combination 
many different background realizations are simulated. The points 
connected with doted lines correspond to the DSFG background 
while the points connected with dashed lines correspond to CMB 
maps. The masses (M200) are reported in Mq. 

We then fit Gaussian functions to the resulting prob- 
ability distributions of the lensed and unlensed back- 
grounds. Figure |3] shows the ratios of the standard de- 
viations of the best Gaussian fits to lensed and unlensed 
distributions {ul/<tij)- This process is repeated for a va- 
riety of lens properties, stepping through mass and red- 
shift. 

4.2. Lensed flux distributions 

Examples of the matched-filtered DSFG and the CMB 
flux distributions are shown in Figure |4) The top (bot- 
tom) panel shows the broadening (narrowing) of the 
DSFG (CMB) flux distribution by lensing, as exp ected 
from the simple arguments presented earlier (j ]2.1 ). The 
high-flux tail of the distribution is caused by occasional 
highly magnified lensing of DSFGs which can substan- 
tially "fill in" the 150 GHz SZ decrement. In all cases, 
lensing conserves the mean of the distributions, keeping 
the average contribution from the DSFG background and 
the CMB unchanged. The lensed DSFG distributions are 
well-fit by the log-logistic function, 



N cx {l3/a){F/a)f^'^[\ + (F/a) 



(11) 



where N is the Poisson expectation value at flux level F, 
a is a scale parameter, /? is a shape parameter. These 
parameters are found to vary as a function of the lens 
properties (cluster mass and redshift). 
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Fig. 4. — Flux distribution of the filtered maps at the location of 
the cluster, in units of matched-filtered map noise. The top panel 
shows the flux distribution of DSFGs while the bottom panel shows 
the distribution of CMB flux in the maps (Af2oo = 9x10^* A/0, 2 = 
0.4). In both panels black curves correspond to unlensed and red 
curves to lensed distributions. The dashed curves show the Gaus- 
sian (bottom) or log-logistic (top lensed) fits to the distributions. 

4.3. Cluster Number Density 

By comparing the lensed and unlensed distributions 
we can estimate the effect of lensing on a cluster survey. 
We compute the number of detected galaxy clusters ex- 
pected in a single-band SPT-like survey, both with and 
without the effects of lensing on the DSFG and CMB 
backgrounds. 

We begin by generating a mass function, a grid of ex- 
pected number densities of galaxy clusters versus mass 
and re dshift, following the prescription of [Tinker et al.| 
( 2012|. Using the signifi cance-mass scaling relation of 



Vanderlinde et al. ( |2010[ ), we convert this to a grid in 
z, where C is the "unbiased significance", defined in 
that work. 

We convert from this grid to an observable space (^— z, 
where ^ is maximum S/N of a detection across filter 
profiles) by convolving in the measurement uncertainty. 
Without lensing of the DSFG background and the CMB, 
this measurement uncertainty has unit width (by con- 
struction, the uncertainty in this space is la). Lensing 
has the effect of slightly increasing the size of this con- 
volution kernel in the case of DSFGs, and decreasing the 
width in the case of CMB lensing. It is important to 
point out that lensing simply changes the variance and 
the shapes of the probability distribution of the back- 
ground emissions while keeping their mean unaltered. 

The expected contributions to noise in matched- 
filtered maps are shown in Figure [5| for both thesingle- 
band SPT cluster survey of Vanderlinde et al. (2010) (the 
exam ple we are using) and t he multi-frequency SP'l' sur- 
2012b[). The unlensed DSFG back- 



vey (Reichardt et al. 



ground contributes 127o of the noise power in the single- 
frequency matched-filtered map, while the CMB con- 
tributes somewhat more, depending on the filter scale. 
At each point in ^— z space, we interpolate an effec- 
tive new level for these contributions based on the ratio 
of lensed-to-unlensed matched-filtered background distri- 
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Fig. 5. — The expected fractional contribution to n oise power in 
the filtered maps for the single-band SPT survey of | Vanderlinde| 
et al. (2010) (top panel) and the multi-frequency SPT survey (e.g., 
Heichardt et al. 2012b, bottom panel), plotted against the core 
size tic ot the p-modci profile used to build the filter. The analysis 
presented in §4 assumes 9c = 0.5'. 

butions (as shown in Figjs]). We convolve this new noise 
level with the mass function, and compare the result to 
the original, effectively isolating the effect of lensing of 
the backgrounds on the observed number count of galaxy 
clusters. The ratio of these number density grids is shown 
by the black curves in Figure |6j for a variety of redshifts 
and significance levels. 
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Fig. 6. — The change in number counts of detected clusters due 
to lensing of the background DSFGs (greater than 1) and CMB 
(less than 1) at various redshifts and detection significances in a 
150GHz single-band SPT cluster survey. 

We find that the lensing of the DSFGs increases their 
effective noise contribution, increasing the Eddington 
bias and preferentially making clusters easier to detect. 
This results in an increase in the number density of de- 
tected clusters. Lensing of the CMB has the opposite 
effect, decreasing the effective noise at cluster locations. 



making it less likely that a cluster would scatter up into 
a fiux-limited sample. 

Finally, we explore the effect of convolving the best- 
fit log-logistic distributions in place of the Gaussian fits 
for the DSFG noise contribution. We interpolate the 
a and /3 parameters (Eq. Ill to each position in the 



^— z grid, calculate the appropriate log-logistic, and con- 
volve it alongside the Gaussian contributions from other 
sources (unlensed CMB, instrumental and atmospheric). 
Taking the ratio of the resulting grid with the original, 
we find a similar set of curves as above, with a slight re- 
duction due to the tail of highly magnified lensed point 
sources, which effectively "fill in" the 150 GHz SZ decre- 
ment. We difference the two sets of lensing-boosted num- 
ber densities to estimate the rate of such strong-lensing 
erasures, and find they occur at roughly the 2% level, as 
shown by the red dashed curves in Figure [61 

We note that while these effects are qualitatively ro- 
bust against the details of the survey being considered, 
the exact magnitude of each will depend on the specific 
properties of the survey, such as frequency coverage, res- 
olution j_j^nd_jioise^^Addiiig_the 90 GHz channel of SPT 



(as in Reichardt et al. 2012b|) reduces the DSFG con 



tribution to the matched-filter noise (see Fig. [s]), such 
that when combined with the CMB contribution the net 
effect of lensing on number counts is near zero. 



5. DISCUSSION AND CONCLUSION 



A s i ngle-b and cluster survey like that of | Vanderlinde| 
et al. (2010) can be affected by strong lensing of either 
■QSFGs or the CMB at the level of roughly 10% in the 
number counts. However, these effects largely offset each 
other, making strong lensing effects a small source of 
bias. While this remains subdominant to Poisson errors 
in the final SPT catalog of ~ 700 clusters (and will be 
further reduced by the multiple frequency bands of the 
full SPT survey), a strong redshift dependence could lead 
to a bias in the inferred growth function or cosmological 
parameter values such as the dark energy equation of 
state. 

The long tail of strongly-lensed background sources 
additionally leads to a ~ 2% decrease in detected clus- 
ter numbers independent of redshift or the original de- 
tectability of clusters. We interpret this to mean that in 
^ 2% of cases, depending only on the geometry of the 
source-lens-observer system, strong lensing of the back- 
ground DSFGs will sufficiently obscure the SZ signal of 
the foreground cluster lens that the cluster would be 
missed. Again, we find this to be a small effect rela- 
tive to the Poisson noise in current catalogs. While this 
is a very dramatic effect when it happens (e.g., SPT- 
CLJ2332-5358 overlaps with a strongly lensed galaxy 
that is obviously obscuring some of the SZ decrement), it 
is extremely rare, and is at most the third-largest strong 
lensing-induced concern (swamped by the increase or de- 
crease in scatter caused by lensing of the DSFGs or the 
CMB). 

These effects will become significant in future surveys 
such as CCAT or SPT-3G, where catalogs containing of 
order 10,000 clusters are expected. They can, however, 
be mitigated with multi-band data, where the contribu- 
tion of the CMB and DSFG background to the final noise 
in the SZ map is reduced. 
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